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In this paper we investigate the cosmological effects of modified gravity with string curvature 
corrections added to Einstein-Hilbert action in the presence of a dynamically evolving scalar field 
coupled to Riemann invariants. The scenario exhibits several features of cosmological interest for 
late universe. It is shown that higher order stringy corrections can lead to a class of dark energy 
models consistent with recent observations. The model can give rise to quintessence without recourse 
to scalar field potential. The detailed treatment of reconstruction program for general scalar-Gauss- 
Bonnet gravity is presented for any given cosmology. The explicit examples of reconstructed scalar 
potentials are given for accelerated (quintessence, cosmological constant or phantom) universe. Fi- 
nally, the relation with modified F(G) gravity is established on classical level and is extended to 
include third order terms on curvature. 

PACS numbers: ll.25.-w, 95.36.+X, 98.80.-k 



I. INTRODUCTION 



One of the most remarkable discoveries of our time is related to the late time acceleration of our universe which is 
supported by observations of high rcdshift type la supernovae treated as standardized candles and, more indirectly, 
by observations of the cosmic microwave background and galaxy clustering. The criticality of universe supported by 
CMB observations fixes the total energy budget of universe. The study of large scale structure reveals that nearly 
30 percent of the total cosmic budget is contributed by dark matter. Then there is a deficit of almost 70 percent; in 
the standard paradigm, the missing component is an exotic form of energy with large negative pressure dubbed dark 
energy Q, 0, y, IH IH 0, 0L The recent observations on baryon oscillations provide yet another independent support to 
dark energy hypothesis Q- 

The dynamics of our universe is described by the Einstein equations in which the contribution of energy content of 
universe is represented by energy momentum tensor appearing on RHS of these equations. The LHS represents pure 
geometry given by the curvature of space-time. Gravitational equations in their original form with energy-momentum 
tensor of normal matter can not lead to acceleration. There are then two ways to obtain accelerated expansion, either 
by supplementing energy-momentum tensor by dark energy component or by modifying the geometry itself. 

Dark energy problem is one of the most important problems of modern cosmology and despite of the number of 
efforts (for a recent review, see 0,0), there is no consistent theory which may successfully describe the late-time 
acceleration of the universe. General Relativity with cosmological constant does not solve the problem because such 
theory is in conflict with radiation/matter domination eras. An alternative approach to dark energy is related to 
modified theory of gravity (for a review, see 0) in which dark energy emerges from the modification of geometry 
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of our universe. In this approach, there appears quite an interesting possibility to mimic dark energy cosmology by 
string theory. It was suggested in refs. [ToL fill [T^ | that dark energy cosmology may result from string-inspired gravity. 
In fact, scalar-Gauss-Bonnet gravity from bosonic or Type II strings was studied in the late universe Eilllll (for 
review of the applications of such theory in the early universe , see 1131 ) . It is also interesting such theories may solve 
the initial singularity problem of the standard big-bang model(see [lj| and refs. therein). Moreover, the easy account 
of next order (third order, Lovelock term) is also possible in this approach (for recent discussion of such gravity, see 

El). 

In this paper we examine string-inspired gravity with third order curvature corrections (scalar-Gauss-Bonnet term 
and scalar-Euler term) and explore the cosmological dynamics of the system attributing special attention to dark 
energy (non-phantom/phantom) solutions. We confront our result with the recent observations. We also outline 
the general program of reconstruction of scalar-Gauss-Bonnet gravity for any a priori given cosmology following the 
method El developed in the scalar-tensor theory. 

The paper is organized as follows. In section two, we consider the cosmological dynamics in the presence of string 
curvature corrections to Einstein-Hilbert action. We analyze cosmological solutions in the FRW background; special 
attention is paid at dark energy which naturally arises in the model thanks to higher order curvature terms induced 
by string corrections. Brief discussion on the comparison of theoretical results with recent observations is included. 
The stability of dark energy solution is investigated in detail. 

Section three is devoted to the study of late-time cosmology for scalar-Gauss-Bonnet gravity motivated by string 
theory but with the arbitrary scalar potentials. It is explicitly shown how such theory (actually, its potentials) may 
be reconstructed for any given cosmology. Several explicit examples of dark energy cosmology with transition from 
deceleration to acceleration and (or) cosmic speed-up (quintessence, phantom or de Sitter) phase or with oscillating 
(currently accelerating) behavior of scale factor are given. The corresponding scalar potentials are reconstructed. It 
is shown how such theory may be transformed to modified Gauss-Bonnet gravity which turns out to be just specific 
parametrization of scalar-Gauss-Bonnet gravity on classical level. Finally, it is shown how to include third order 
curvature terms in the above construction. Summary and outlook are given in the last section. 



II. DARK ENERGY FROM HIGHER ORDER STRING CURVATURE CORRECTIONS 



In this section we shall consider higher order curvature corrections to Einstein-Hilbert action. To avoid technical 
complications we restrict the discussion up to third order Riemann invariants coupled to a dynamical field </>. The 
cosmological dynamics of the system will be developed in detail and general features of the solutions will be discussed. 
It is really interesting that the model can account for recent observations on dark energy. 



A. General action 



We begin from the following action 



S 



d X\J—g 



— 2 - \w{(t>)g^d^d v ^ - V(cj>) + C c + C 



2 k 



(1) 



where is a scalar field which, in particular case, could be a dilaton. C m is the Lagrangian of perfect fluid with 
energy density p m and pressure p m . Note that scalar potential coupled to curvature (non- minimal coupling) |l7j does 
not appear in string-inspired gravity in the frame under consideration. 
The quantum corrections are encoded in the term 

£ c = awo4 1} +6(0)4 2) (2) 

where £i(0) and ^(^) are the couplings of the field <fi with higher curvature invariants. C c and ci^ are given by 



C^=E 3 +R^Rf p R^ 
The third order Euler density E$ is proportional to 

Ei oc f > lv P< TTr > f ,,,,,,/? p' v ' R p' & ' U r '''' 

±-13 u< c e^< pi ' a i ' T ' 'rj' -il^j, ll pa lx, Tr\ 



(3) 
(4) 



(5) 
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Since there does not exist t^ l/paTri if the space time dimension D is less than 6; £3 should vanish when D < 6, 
especially in four dimensions. By using 

e^V.'pV'rV = 6/5/6/5/6/6* ± (permutations) , (6) 
we can rewrite the expression ||SJ| as 

B 3 a 8 (R 3 - URR^R^ + 3RR^ pcT R^ pa + WR^RfR^ + 24R" R p a R u f p 

-24r;r^rjp + 2R^R P J v R T r - WKT) • ( 7 ) 

We should note in the r.h.s. of ©, there appears 6! = 720 terms, which correspond to the sum of the absolute values 
of the coefficients in each term in the RHS of (Q) 

8(1 + 12 + 3 + 16 + 24 + 24 + 2 + 8) = 720 . (8) 

In what follows we shall be interested in the cosmological applications of modified equations of motion and thus 
assume a flat Friedmann-Robertson- Walker (FRW) metric 

d 

ds 2 = -N 2 (t)dt 2 + a 2 (t) J2( dxl f, ( 9 ) 

i=l 

where N(t) is the lapse function. With the metric ©, the Riemann invariants read 

4" = 24* 2 (*±^ - £f) , 4 2) = ^ + I») _ Z^W (10 ) 

where I = H + H 2 and iJ = d/a. It is straightforward though cumbersome to verify explicitly that third order Eulcr 
density E 3 is identically zero in the FRW background. The non- vanishing contribution in Ea. i|10f> comes from the 
second term in Q. To enforce the check in a particularly case, we consider D dimensional de-Sitter space, where 
Riemann curvature is given by 

R^r = Ho (<yv - w) ■ (11) 

Here Hq is a constant corresponding to the Hubble rate. In the de-Sitter background we have 

E 3 cx D{D - 1){D - 2){D - Z){D - 4){D - 5) , (12) 

which is obviously zero in case of D < 6. For simplicity we shall limit the discussion to a homogeneous scalar field 
4>{t). Then the spatial volume can be integrated out from the measure in equation l|T]l. which we rewrite as 

S = J dtNa 3 A. + C c + + C m . (13) 

where = —^uj((f>)(V(j)) 2 — V (</>). Varying the action Ijl3(l with respect to the lapse function iV we obtain jljj 

— 5- = Pm + + Pc (14) 

where 

P<{, = + ^(0) (15) 

In Eq.lO, thc 

energy density p c is induced by quantum corrections and is given by the following expression 



, dC c d dC c dC c 
Pc= 13 %v" + ^dA^ diV ' 



(16) 



JV=1 



It would be convenient to rewrite p c as 

Pc = Hi^)p£ ) +U<f>)P { c ) (17) 
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Using Eas. (|10f) & l|10(l we obtain the expressions of p^ and p^ 

pf> = -24ff 3 Si (18) 
p (2) = -72HI 2 E 2 - 72 (ill 2 + 2IHl\ - 216H 2 ! 2 + 120 (H 6 + I 3 ) (19) 

where Si = £i/£i and S2 = £2/62 It is interesting to note that the contribution of Gauss-Bonnet term (described 

by Eq.^niO cancels in equations of motion for fixed (f> as it should be; it contributes for dynamically evolving scalar 

field only. In case of the third order curvature corrections, the Euler density is identically zero and hence it does not 

(2) 

contribute to the equation of motion in general. Secondly, C c contributes for fixed field as well as for dynamically 
evolving <p. It contains corrections of third order in curvature beyond the Euler density. 

We should note that such higher-derivative terms in string-inspired gravity may lead to ghosts and related instabil- 
ities (for recent discussion in scalar-Gauss-Bonnet gravity, see 18] ). However, the ghost spectrum of such (quantum 
) gravity (for the review, see^||) is more relevant at the early universe where curvature is strong, but less relevant at 
late universe. Moreover, in accordance with modified gravity approach, the emerging theory is purely classical, effec- 
tive theory which comes from some unknown gravity which has different faces at different epochs. (Actually, it could 
be that our universe currently enters to instable phase). For instance, in near future the currently sub- leading terms 
may dominate in the modified gravity action which then has totally different form! Hence, this is that (unknown) 
gravity, and not its classical limit given by Eq.Q relevant during specific epoch, whose spectrum should be studied. 
The point is best illustrated by the example of Fermi theory of weak interactions whose quantization runs into well 
known problems. Finally, on the phenomcnological grounds, it is really interesting to include higher order terms. At 
present the situation is remarkably tolerant in cosmology, many exotic constructions attract attention provided they 
can lead to a viable model of dark energy. 

The equation of motion for the field <fi reads from l|13fl 

uj(4> + 3H<j>) +V- - <e 2 4 2) + u<j> - Lj'-j- = (20) 

In addition we have standard continuity equation for the barotropic background fluid with energy density p m and 
pressure p rn 

p m + 3H(p m + Pm )=0 (21) 

Equations (|14|) . H20|) . and (|21|l are the basic equations for our system under consideration. 
Let us note that in the string theory context with the dilaton field (f> we have 

V{cj)) =0, & = cia'e 2 ^ , £ 2 = c 2 a' 2 e 4< ^ (22) 

where (ci, 02) — (0, 0, 1/8), (1/8, 0, 1/8) for Type II, Heterotic, and Bosonic strings, respectively. 



B. Fixed field case: general features of solutions. 

We now look for de-Sitter solutions in case of <fi = constant and p m = 0. In this case the modified Hubble Eas. i|14|) 
gives rise to de-Sitter solution 

/ 1 \i/4 

3 = 24£ 2 # 4 or H= ^— J (23) 

where £2 = ^exp(— 4^/0o) for type II and Bosonic strings. Normalizing £2 to one, we find that H = 0.6 (we have 
set k 2 — 1 for convenience). Below we shall discuss the stability of the solution. There exists no de-Sitter solution 
for Heterotic case. Actually, de-Sitter solutions were investigated in the similar background in Ref.^3] where higher 
order curvature corrections up to order four were included. Since, here we confine ourselves up to the third order and 
the fourth order terms are excluded from the expression of p c ; these terms come with different signs. Thus it becomes 
necessary to check whether or not the stability property of de-Sitter solutions is preserved order by order. 

We further note that the modified Hubble Eas. H14|l admits the following solution in the high curvature regime in 
presence of the barotropic fluid with equation of state parameter w 



a(t)=a t h °, or a(t) = a (t s - t) hl 



(24) 
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where 

(25) 



a 



1 + w ' 



(72(-Mo + 2J„Mb) + 216^/ 2 - I20(ftg + 7 3 )) 



(26) 



We have used p rn — p$a 3 ( 1+tu ) for the background matter density and Iq = ho(h — 1), I = — 2/i (^o — !)• F° r the 
effective equation of state dictated by the modified Hubble Eas. <|14fl we have 



2 i? 1 + w 

3H 2 " - ~ 1 + ~ 

It is interesting to note that effective EoS parameter l|27|l may correspond to inflationary solution in the presence 
of background fluid (radiation/matter). In the low curvature regime or at late times w c s — w. In the presence of 
phantom matter, the effective EoS being less than — 1 is typical for Big Rip singularity. It is really not surprising that 
we have inflationary solution at early epochs in the presence of higher order curvature correction to Einstein Hilbert 
action; an early example of this phenomenon is provided by R 2 -gravity. 

C. Autonomous form of equations of motion 

Let us now cast the equations of motion in the autonomous form. Introducing the following notation (k 2 = 1) 

x = H, y = H, u= 4>, v = <j>, z = p m (28) 
We shall assume u>(4>) = v = const, we obtain the system of equations 
x = y, 

b 



)vv 2 - 24£iSi:e 3 + \-72xI 2 E 2 - 72(yl 2 + 4Iyx 2 ) - 216x 2 I 2 + I20(x 6 + I 3 )} - 3x 2 

V = 



U4I(x,y)&x lMI&x' 

u = v , v = . z — — 6x(l + w)z (29) 

We shall be first interested in the case of fixed field for which we have (assuming v = 1) 

x = y, 

-72(yl(x, y) 2 + 41 (x, y)yx 2 ) + 2\&x 2 I 2 {x 1 y) + 120{x 6 + I 3 (x, y)) - 3a; 2 + z] 
JJ ~ U4I(x,y)x ' 

z = -3x(l + w)z (30) 

where 

I{x,y)=x 2 +y, = 24x 2 (y + x 2 ), = 24(x 6 + I 3 {x 7 y) (31) 

In the presence of perfect fluid, the de-Sitter fixed point is characterized by 

s c = 0.71, y c = Q, z c = Q (32) 
Perturbing the system around the critical point and keeping the linear terms we obtain 
Sx = Sy , 

/21 2 10 1 \ _ (2 5 1 \ 1 . 

by = ^rx c + — + -—5 )Sx+[-x e + — s + —3 Sy + —-jSz , 



^3 c 3x c 48x 2 J \3 6x1 ^ 3 J y U4x 3 
Sx = — 3x c (l + w)Sz (33) 
Stability of the fixed points depends upon the nature of eigenvalues of perturbation matrix 

1 



Ai,2 = ~ ( «22 ± V 4a 2 i + a\ 2 ) , A 3 = a 33 = -3x c (l + w) (34) 



For the fixed point given by Ij32(l . Ai is positive where as A2 is negative making the de-Sitter solution an unstable 

node. In 
FIG.ITJ. 



node. In fact, Ai remains positive for any x c > thereby making the conclusion independent of the choice of £2°' (see 
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FIG. 1: Plot of the first eigenvalue Ai versus the critical point x c . The eigenvalue remains positive if the critical point is varied 
from zero to larger values. 




FIG. 2: Plot of Y = (^S~) x 10 5 (black line) and Y = (^£-) (gray line) versus X = 0g. corresponding to ho = 40 or 

wde = —0.95 (y is assumed to be one). The common region corresponding to positive values of the couplings gives possible 
models of dark energy induced by higher order curvature corrections. 



(i) 



(2) 



D. Dynamically evolving field <f> and dark energy solutions 

In what follows we shall be interested in looking for an exact solution of equations of motion (|14fl and (|2U[1 which 
of interest to us from the point of view of dark energy in absence of the background fluid. In this case let us look for 
the following solution 

H = — , <j> = </>o In — - (when h > 0) , H = ^° , 4> — <t>o l n — — - (when ho < 0) . (35) 
t t\ t s — t t\ 

Substituting (|35|1 in evolution Eqs. (|14|l and H2U[) yields (we again set k 2 = 1) 



1/(1 - 3ho)4 + -^hh 3 (h Q - 1) + -g-(h% + / 3 ) = o , 

-3^ + - + = (36) 



where 



J = ^ (-288^0 - 72(-/i / 2 + 2/oAoio) - 216/i /o + 12 0(^o + *o)) . 

/o = fto(fto - 1), /o = -2/i (/io - 1) (37) 



7 



0.1 r- 

0.09 " 




FIG. 3: Plot of Y = x 10 5 (gray line) and Y = (^M (black line) versus X = 0§. corresponding to ft = -33.33 or 

wde = —1.06. The region bounded by 6 < X < 31.5 corresponds to possible phantom dark energy models. 



Using Eqs.lOSJ), we express the couplings through ho and <po 



(0) 



96# 



3^ 



^ + ^ _ 1)02 



1 

"o 



J(h )(h -l) + h 3 (h 3 + (h -l)^) 

ih 2 a - v -f + v(3h - 1)$ ) J(h ) 



-m + + 



J(M(^o-l) + ^(^ + (/io-l) 3 ) 



Let us note that the string couplings (£i(0) 



^,6(0) = e 



(0) c 



with m 



(38) 
(39) 

4) are generic to solution 



described by (|33|) : for other couplings such a solution does not exist. We also note that Eas. l|55|l & l|3^|l reduce to 
the earlier obtained results in Ref . 01 ( see Refs. [l3L I20L l2l| on the related theme) where similar investigations were 
carried out confining to only second order curvature invariants in the action Q . 

There are several free parameters in the problem. In order to extract important information from Eqs. (|38|1 and 
(|39|l . we proceed in the following manner. We fix ho corresponding to the observed value of dark energy equation of 



(0) 



and £<°> 



leading to allowed values of 



that for a given value of h from the 
Our model incorporates higher order curvature corrections 
This gives rise to comfortable choice of the equation of state 



state parameter wde and impose the positivity condition on the couplings i 
the parameter <pQ. In the absence of coupling $2(0)7 it wa s shown in Rcf. 1 
allowed interval, the parameter O takes a fixed value, 
allowing a one parameter flexibility in the values of </>o- 
consistent with observations. 

The three years WMAP data is analyzed in Ref.[z|. which shows that the combined analysis of WMAP with 
supernova Legacy survey (SNLS) constrains the dark energy equation of state wde pushing it to wards the cosmological 
constant. The marginalized best fit values of the equation of state parameter at 68% confidence level are given by 
— 1.14 < wde < —0.93. In case of a prior that universe is flat, the combined data gives —1.06 < wde < —0.90. Our 
model can easily accommodate these values of wde- For instance, in case of non-phantom (standard dark energy) 
we find a one parameter family of dark energy models with ho ~ 40 {wde = —0.98) corresponding to 0g > 41. Like 
wise, in case of phantom dark energy, we find that for ho — —33 (wde = —1.02), the viable range of the parameter 
(j>Q is given by 6 < 0p < 31.5. These values are consistent with recent WMAP release and SNLS findings. 

We should mention that the observations quoted above do not incorporate the dark energy perturbations which 
might severely constrain the phantom dark energy cosmologies. The combined data (CMB+LSS+SNLS) then forces 
the dark energy equation of state to vary as, —1.001 < wde < —0.875 22J. Our model can easily incorporate 
these numerical values of wde by constraining ho and </>g similar to the case of non-clustering dark energy. A word 
of caution, the evolution of dark energy perturbation across the phantom divide needs additional assumptions; a 
complete analysis should take into account the non-adiabatic perturbations which makes dark energy gravitationally 
stable 0. 
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Dark energy 


ho 


4>l 


WDE 


Observational constraint on wde 


Constraint on wde with flatness prior 


Non-phantom 
Phantom 


40 
-33.33 


6 < 4>\ < 31.5 


-0.98 
-1.02 


-1 06+ 013 


u - 3 ' -0.09 



TABLE I: Observational constraints on (non-clustering) dark energy equation of state wde dictated by the combined analysis 
of WMAP+SNLS data[22l and the numerical values of model parameters consistent with the observations. 



E. Stability of dark energy solution 



In what follows we shall examine the stability of the dark energy solution 1|35|) induced by purely stringy corrections. 
In general the analytical treatment becomes intractable; simplification, however, occurs in the limit of large ho 
corresponding to w e f / ~ — 1. 

Let us consider the following situation of interest to us 



Prn = 0, W = V = 1, 6(0) = CfV^ , 6(0) = d ^' 



(0) 40/00 



(40) 



In order to investigate the stability around the dark energy solution defined by l|35|) . we need a convenient set of 
variables to cast the evolution equations into autonomous form. We now define the variables which are suited to our 
problem. 



X 



H 



y=[H + H') 6(0) , Z = H%(4>) , 



did 
dN = Hdi 



With this choice, the evolution equations acquire the autonomous form 



dX 
dN 



dy_ 

dN 

dZ 
dN 



— - - l 'Kl 0) ( - -:ij>> 



X 48 



y 



1/2 



96 e 



(0) 



mi / \ V 2 



ix 2 2 2f\ 0) y ( y 

— _ zx -2y + SLJL I JL 

24k 2 144 30o 



1/2 



X \Z + 26 



(o) 



AO) 

A 2 , 



3Z 

1/2 



AO) 

S2 , 



b Z 



5 A0) 



AO) 

-S2 , 



(41) 



(42) 



We have used the field equation (|2(J|I and Ea. (|29|l for H in deriving the above autonomous form of equations. For our 
solution given by (|35f) . we have 



X — Xq 



00 



h 2 f (0) 

Z = A) = 



t 2 
' l 



y = y = 



(-ho + h 2 ) 2 ^ 



t 4 
1 1 



(43) 



It can be checked that (X ,y , Z ) is a fixed point of ijUJl. We then consider small perturbations around (|43|l or 
equivalently around the original solution (|35|l 

x = x Q + sx , y = y a + sy , z = z + sz . (44) 

Substituting (|44|l in (|42|l and retaining the linear terms in perturbations, we find 



fsx\ 




f SX 


sy 


= M 


sy 


[sz) 




{ sz 



d 

dN 



Here M is 3 x 3-matrix perturbation matrix whose components are given by 

„ , -l + h 



(45) 



Mil 



h 
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Mi 



<M? 24£< 0) t? , 144(-1 + V) 



M 



13 



M 2 1 

M 22 = -2- 



M 23 

M31 
M 32 



2^o 0<4 O) (-/io + ftg) Mo 
<M?(-1 + M t 192^ 0) hg 96(-l + fe ) 3 d 0) 

1 gfe^ 
72 3</> i? ' 
(-1 + M 

2 2(-l + feo) 3 gf 10fe^ 0) 



AO). 2 
si c i 



^ 0) (-fto + 

M 33 = -2 + A . (46) 

Stability of the fixed point(s) depends upon the nature eigenvalues of the perturbation matrix M. If there is an 
eigenvalue whose real part is positive, the system becomes unstable. Here for simplicity, we only consider the case of 
h — > ±00, which corresponds to the limit of w e g ~ — 1. In the case, we find 



sT 1 1 



t\ 40h§ ' t\ 32h 

and the eigenvalue equation is given by 



(47) 



= F(A) = -A 3 -6A 2 --^_A-^2_ . (48) 
v ' 4O0 O 4O0o 

The values of A satisfying F[X) = give eigenvalues of M. The solutions of (|48(l is given by 



A = A± = ± 



+O(\ho\) , A = A = — 7 + (|^or 1/2 ) • (49) 



40 



When ho < 0, the mode corresponding to A+ (A_) becomes stable (unstable). Since A± are pure imaginary when 
ho > 0, the corresponding modes become stable in this case. On the other hand, the mode corresponding to Ao is 
always stable. Thus, the non-phantom dark energy solution (|35(l induced by string corrections to Einstein gravity 
is stable. Such a solution exists in presence of a dynamically evolving field <j) with V(<f>) — coupled to Riemann 
invariants with couplings dictated by string theory. Dark energy can be realized in a variety of scalar field models by 
appropriately choosing the field potential. It is really interesting that we can obtain dark energy solution in string 
model without recourse to a scalar field potential. 

Let us compare the results with those obtained in [Hj], where £ 2 = but V{4>) 7^ 0. The dark energy solution 
studied in Ref.|lOj|. was shown to be stable when ho > but unstable for ho > 0. The present investigations include 
£ 2 and V = which makes our model different from Ref. 10]; it is therefore not surprising that our results differ from 
Ref . [Tol] . Since ho > corresponds to the quintessence phase and ft, < to the phantom, the solution in the model 
(with £0 and V — 0) is stable in the quintessence phase but unstable in the phantom phase. We should notice that 
the approximation we used to check the stability works fine for any generic value of ho- For instance, 5 < ho < —667 
which corresponds to the variation of wde in case the dark energy perturbations are taken into account. We also 
carried out numerical verification of our results. 



III. THE LATE-TIME COSMOLOGY IN SCALAR-GAUSS-BONNET GRAVITY 



A number of scalar field models have recently been investigated in connection with dark energy (see Ref.Q for 
details). The cosmological viability of these constructs depends upon how well the Hubble parameter predicted by 
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them compares with observations. One could also follow the reverse route and construct the Lagrangian using the 
observational input; such a scheme might help in the search of best fit models of dark energy In what follows we 
shall describe how the reconstruction program is implemented in presence of higher order string curvature corrections. 



A. The reconstruction of scalar-Gauss-Bonnet gravity 



In this section it will be shown how scalar-Gauss-Bonnet gravity may be reconstructed for any requested cosmology 
using the method |ljj developed in the scalar-tensor theory. We limit here only by Gauss-Bonnet term order (by 
technical reasons) but there is no principal problem to include higher order terms studied in previous section. It is 
interesting that the principal possibility appears to reconstruct the scalar-Gauss-Bonnet gravity for any (quintessence, 
cosmological constant or phantom) dark energy universe. The last possibility seems to be quite attractive due to 
the fact 01 j that the phantom universe could be realized in the scalar-Gauss-Bonnet gravity without introducing 
ghost scalar field. In this section, we show that in scalar-Gauss-Bonnet gravity, any cosmology, including phantom 
cosmology, could be realized by properly choosing the potential and the coupling to the Gauss-Bonnet invariant with 
the canonical scalar. 

The starting action is 



S = d X\f—g 



(50) 



Here G is the Gauss-Bonnet invariant G = Cc © and the scalar field <j> is canonical in 1)500. As in previous section, 
it is natural to assume the FRW universe 10 with N(t) = 1 and the scalar field <j> only depending on t. The FRW 
equations look like: 

= --H 2 + -4? + V{<f>) + 2AH 3 >' , 51 

k z 2 at 

0=\(2H + 3HA + V - Vifi) - 8i^(f)) - WHH^m - 16ff 3<W)) . (52) 
k V J 2 at 2 at at 

and scalar field equation: 

= 4> + 3H<t> + V'((f>) + &(<f>)G . (53) 

Combining 1)51(1 and 1)52(1. one gets 



Ea.l(5i|l can be solved with respect to £i((j>(i)) as 



Combining ((51(1 and l)55|). the scalar potential V((f>(t)) is: 



3 TTI 2 1 ; /iN 2 n /|NTW|N f dt l ( 2 tt,, \ , I,, ^ 2 



VW)) = ^H(ty - -4>{t) -3a(t)H(t)J -^-l—H(ti) + <l>(h) j . (56) 

We now identify t with f{4>) and H with g'(t) where / and g are some unknown functions in analogy with Ref.|l6| 
since we know this leads to the solution of the FRW equations subject to existence of such functions. Then we consider 
the model where V((f>) and £i(0) may be expressed in terms of two functions / and g as 

V&) = ^9'(fm 2 -^,-3g'(m)e^ f Vw-^ x + , 

ft( *> = I / / *"* )f11 *° b 9 " <««> + ?w) ' (57) 

By choosing V(<p) and £i(0) as 1)57)1 . we can easily find the following solution for Eas. ()51l) and 1)52(1: 

4>=f- 1 {t) (t = /(0)) , o = o e"W (H = g'(t)) . (58) 
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We can straightforwardly check the solution l|58|l satisfies the field equation (|53|) . 

Hence, any cosmology expressed as H — g{(f)) in the model l|5U|l with (|57|) can be realized, including the model 
exhibiting the transition from non-phantom phase to phantom phase without introducing the scalar field with wrong 
sign kinetic term. 

In the Einstein gravity, the FRW equations are given by 

= -^H 2 + p, 0= \ (2H + 3H 2 ) +p . (59) 

Here p and p are total energy density and pressure in the universe. By comparing (|59() with l|58[l wc find the effective 
energy density p and the pressure p are given 

p = , P = \g"{t) ■ (eo) 

Since t = g' 1 ((k) y/ p/Sj , we obtain the following effective equation of the state (EoS): 

P = -p-\g"(g'- 1 (^)) , (6i) 



k 2 \ V V 3. 

which contains all the cases where the EoS is given by p — w(p)p. Furthermore, since g' 1 could NOT be always a 
single- valued function, Ea. (|61[l contains more general EoS given by 

= F(p,p) . (62) 

This shows the equivalence between scalar-tensor and ideal fluid descriptions. 

Let us come back now to scalar-Gauss-Bonnet gravity. It is not difficult to extend the above formulation to include 
matter with constant EoS parameter w m = p m j p m . Here p m and p m are energy density and pressure of the matter. 
Then, instead of lfST|) and (|5*2~)l the FRW equations are 



= -\h 2 + U 2 + V{4>) + Pm + 24i/ 3 ^M , (63) 
n z 2 at 



« = \(2H + 3H 2 ) + ^ - V{4>) +Pm - S^^PI - 16HH^1 - Iftfl-^pi . (64) 
k V / 2 ' at at at 

The energy conservation law 

p m + 3H (p m + p m ) = , (65) 

gives 

Pm = PmOCT a(1+ ™ m) , (66) 

with a constant p m o- Instead of 1)57(1 . if we consider the model with 

+ (l + Wm ). 9o e- 3 ( 1 +-")^^))) , 
= d ^ 7 ^/'(0 2 )e^^»^ 5 "(/(0 2 )) 

+ 2^ + (l + - m ) 5 oe- 3 ( 1+ ^^^))) , 
we re-obtain the solution l|58|) even if the matter is included. However, a constant ao is given by 



(67) 



a = - ■ (68) 
Po 
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One can consider some explicit examples |l6[ : 



« = /(0) = -T' ff(t) = holnT 7 , 

00 Is - t 



which gives 



H = h 



t U 



H 



h t s (2t - t a ) 
t 2 (t s - t) 2 



(69) 



(70) 



Then the universe is in non-phantom phase when t < t s /2 and in phantom phase when t > t s /2. There is also a 
Big Rip singularity at t = t s . Especially in case w m = (that is, matter is dust) and ho = 2, we reconstruct the 
scalar-Gauss-Bonnet gravity with following potentials: 



v{<j>) = 



a wo 



Here 



60o0 




40^0 S 




-( 


K 2 0(0 S - 


- 0) - 


{4>.-<f> 


f 


K 2 I 


+ 50 (_ 


# 401 


280^ 


1405 


70 7 30 6 


50 5 




4 


1 


f 4 Z0202 






2 4 


320§02 


K 2 \ 6 






4 


+ 140^0 


35040 2 
3 


28030 3 

3 


- 70^0 4 



30 3 2 

70a!. 4 

+ 



20 s In -f- 



303 

- 20 s 



5 S 25 



5 0, 25 



In — 



</> 6 , 5o / 



Ci 



6 0g V 14 '/' 2 



401 2806 ^ 
30 5 0s 



C2 



(71) 



)£ s and ci, C2 are constants of the integration. 



Another example, without matter (go = 0), is [2c 



COS C/Q 



sinwi , f-\t) 



UOt 

b sm— . 



(72) 



Here ho, 9o, to, and 0o are constants. This leads to reconstruction of scalar-Gauss-Bonnet gravity with 



V{<t>) 

Then from Eq.J72J we find 



3/iQ 
*-2 



1 + COS ( 



2cos#o ,2^ ^o^ 2 



QJ<Po 
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1 + cos ( 



2\ 1/2 



2 cos i 



^2 
-<Pl 



H = ho (1 + cos 6>o cos wi) > , £f = —hou> cos #o sinwi , 



(73) 



(74) 



Then the Hubble rate H is oscillating but since H is positive, the universe continues to expand and if hococos9o > 0, 
the universe is in non-phantom (phantom) phase when 2nn < tot < (2n + 1) tt ((2n — 1) tt < ujt < 2nir) with integer 
ii. Thus, the oscillating late-time cosmology in string-inspired gravity may be easily constructed. 
One more example is [23( 



TJ 

g(t) = Hot lncoshi?o^ 

Ho 

Here we assume H > Hi > 0. Since 

H = g'(t) = Ho-H 1 tanhHot , H = g"{t) = 



Ho Hi 
cosh 2 H t 



< 



(75) 



(76) 



when t — ► ±oo, the universe becomes asymptotically deSitter space, where H becomes a constant H — > ifg T -ffi and 
therefore the universe is accelerating. When t = 0, we find 



- = H + H 2 = -HxH + H 2 < , 
a 



(77) 
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therefore the universe is decelerating. Then the universe is accelerating at first, turns to be decelerating, and after that 
universe becomes accelerating again. As H is always negative, the universe is in non-phantom phase. Furthermore 
with the choice 



0, i = /(</>) 



1 ( kH 
— tan — = — . 
Hq \2y/2H 1 



(78) 



we find the corresponding scalar-Gauss-Bonnet gravity 

Hi 



v(4>) 



2 

— (H - Hi tanh tp) 



^ {H - Hi tanh^) (l + e 2v 



Hn 



,9o 
2Ho 



dip' 



\/2k 2 cosh 2 

2p- 

V 



In (1 + c 2 ^) + 



(i?o - ffi tanh</) z 



Here 



2<p' - In 1 



9? = tan 



kHo 



6(1 + e 2 ^) 
5 



6(1 + e 2 vy 
5 



6(1 + e 2 v) c 
2 



6(l + e 2 ^') 6(l + e 2 ^') 6 (1 + e 2 ^')' 



(79) 



(80) 



Although it is difficult to give the explicit forms of V{<j>) and £i (</>), we may also consider the following example 



/(/) = //„ ( l - 



h + t 2 ^ ti + t 2 ^ 2 



6 2 

Here ho, ti, t 2 are constants. Hence, Hubble rate is 

H(t) = h (j- t -±±^t 2 + t 1 t 2 t 



(3ti > t 2 > h > , ho > 0) 



H(t)=h (t-h){t-t 2 ) 



(82) 



Since H > when i > and H < when i < 0, the radius of the universe a = aoe 9 ^ has a minimum when 
t = 0. From the expression of ij in lj%2l . the universe is in phantom phase (h > O^j when t < t\ or t > t 2 , and 

in non-phantom phase (^H < 0^ when ti < t < t 2 (for other string-inspired models with similar cosmology, see for 

instance Then we may identify the period < t < ti could correspond to the inflation and the period t > t 2 

to the present acceleration of the universe (this is similar in spirit to unification of the inflation with the acceleration 
suggested in other class of modified gravities in ref. 25]). If we define effective EoS parameter w e s as 



P , 2H 
p SH Z 



(83) 



we find w c s — > — 1 in the limit t — > +oo. Although it is difficult to find the explicit forms of V(<fi) and £i (</>), one 
might give the rough forms by using the numerical calculations. From the expression of V{4>) and £i(</>) in (|57fl . if 
f(<p) is properly given, say as t — f((f>) — 4>/4>o with constant <po, there cannot happen any singularity in V(cj>) and 
£i(4>) even if t — t\ or t = ti, which corresponds to the transition between phantom and non-phantom phases. Then 
the model (I5U[) could exhibit the smooth transition between phantom and non-phantom phases. 
The next example is 



j(t) = h In 



Since 



to 



H 



t = f(<f>)=t et 



ho 
t 



(84) 



(85) 



we have a constant effective EoS parameter: 



w cB = - 



(86) 
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Eqs.JEIJl give 

v(4>) 
CM 



f3h 2 (l-h ) , 4>l 



(ho + l)t% V 



16/ig (fo + i) 



K 2 



1 - 5/iq) e 



■i± , 3/l (l+W m )50 _3(l + "-m)HQ0 

*o -| : e *o 

(4 + 3w m ) ho - 1 



+ 0o e*o 



+-{3(l+w m )/i -4} 1 {(4 + 3w m ) /i - 1} 1 (1+tWm) ffotoe" 



{3(l + mm )h -4}, 
00 



(87) 



Thus, there appear exponential functions, which are typical in string-inspired gravity. 

As clear from Q85|l. if /iq > 1, the universe is in quintessence phase, which corresponds to —1/3 < w e g < — 1 in 184fl . 
If ho < 0, the universe is in phantom phase with w e s < — 1. In phantom phase, we choose to to be negative and our 
universe corresponds to negative t, or if we shift the time coordinate t as t — > t — t s , with a constant t s , t should be 
less than t s . 

The model ^3 corresponds to g a = in JS7J|. In the notations of ref.0]' to = ii, ^(</>) = Vbe~*o , and /(^) = 

20 

/oe^o = — £i(0). Then from the expression l|87|) . one gets 



To 



1 / 3/i§(l-ft ) $ , 



/o 



2ft 

16ft§ (/i + 1) V 



which is identical (after replacing to with ti) with (16) in |l(J |. 

Thus, we demonstrated that arbitrary late-time cosmology (from specific quintessence or phantom to oscillating 
cosmology) may be produced by scalar-Gauss-Bonnet gravity with scalar potentials defined by such cosmology. The 
reconstruction of string-inspired gravity may be always done. Moreover, one can extend this formulation to include 
the higher order terms in low-energy string effective action. 



B. The relation with modified Gauss-Bonnet gravity 

In this section we show that scalar-Gauss-Bonnet gravity may be transformed to another form of modified Gauss- 
Bonnet gravity where no scalars present. In addition, the formulation may be extended to include higher order terms 
too. Starting from l|5U|l . one may redefine the scalar field <f> by <j> — etp. The action takes the following form 



S = I d 4 a;v '—g 



Here 



V(<p) = V(e<p), fcM. 
If a proper limit of e — ► exists, the action (|89() reduces to 



S = I d^x^j—g 



*-V(sp)-U<P)G 



Then (p is an auxiliary field. By the variation of 99, we find 

= V"(¥>)-fi(¥>)G 

which may be solved with respect to ip as 

<p = $(G) . 

Substituting (JSIJ into the action the F(G)-gravity follows pt| : 



S 



d^x^f—g 



R 
2^ 



F(G) 



, F(G) = V-($(G))-|i(*(G0)G 



(89) 

(90) 

(91) 

(92) 
(93) 

(94) 
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For example, in case of Q84J1. in e — > limit after redefining <j> = ecp and <fio = epo, V(<p) and £i(4>) reduce to 
3/igOo-l) . 3h (l + w m )g _3(i+^ m) h O 



(/i + 1) t 2 ^ 



e + 



(4 + 3w m ) ho - 1 



+ £ {3 (1 + w m ) ho - 4} 1 {(4 + 3u> m )/»o - 1} 1 (1 + w m ) gotfc ^^P' ~ .(95) 



{3(l + mnl )h -4H 



8/i (/to + 1) k 2 8 
The solution corresponding to (|84ll is 



g(t) = ho In ^- , 99 = 99 ln— . 

to to 



If we further consider the case <?o = 0, Eq. l|92|l gives 



e 



/J 



24/ig (/i - 1) 



G 



(96) 



(97) 



Eg. (197(1 could have a meaning only when ho > 1 or ho < if G is positive. In this situation 

1 



F(G) = AoG 1 ' 2 , A 



3 (h - l )h 



2(1 + /i )k 2 V 2 

The above model has been discussed in [2(j . Actually, in (2(| the following type of the action has been considered: 



S = J d 4 x^y^R + F(G) 



(99) 



In case that F(G) is given by JHHJ, m terms of [2q|. A = fo- Hence, A Q (jHEl coincides with the Eq.(26) of |2q . 

As a further generalization, we may also consider the string-inspired theory of second section where next order term 
is coupled with scalar field: 



(100) 



As in JH3>, 

we may redefine the scalar field <pby <j) = ep. If a proper limit of e — > exists, the action (|100|l reduces to 

i? 



2 k 



Here 



6 = hW 2 (ev?) ■ 

e— »0 

Then ip could be regarded as an auxiliary field and one gets 

= V'(p)-£[(p)G + &£W 

which may be solved with respect to p as 

p = ^{g,c { < 2) 

Substituting (|104|) into the action l|101|) . we obtain F(G, Cc )-gravity theory 
R 



(101) 



(102) 



(103) 



(104) 



S = I d^xyj—g 



, ^(g,4 2) ) 



2k 



In case of the string-inspired gravity: 



, F(G) = V ($(C?,4 2 ))) ($(G,4 2) )) G + & (<f(G,4 2) )) 4 2) ■ (105) 



2<j> 2<j> A£ 

V = V e ^ , £ 1= £ e<»o , ^2 = %e 7 o 



(106) 
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Here <fio, Vq, £o, and rjo are constants. We may consider the limit of e — > after redefining <fi = eip and 4>q — e</?o- 
Thus, Ea. fTU3|) gives 



e^ = e(G,4 2 ))^-M_ + r(G,4 2) ) 



Here 



and 



Y(G,£? ) ) = y++: 



y± 



v a 



(2) 



± 



y+ea" +y_e! 



1/3 



V 



(2) 



6rj £ 



(2) 



(107) 



(108) 



(109) 



(2) 

Hence, the action of the corresponding F(G,£ C )-theory is 



S = I d 4 Xy/—g 



F(G,4 2) ) 



V 



fr 2 -F{G,c { ?) 

- £ e (g, 4 2 >) g + ??0 e (g, 4 2) )' £ 



(2) 



(110) 



e(G,d 2) ) 

Instead of (|50|l . one may consider the model with one more scalar field x coupled with the Gauss-Bonnet invariant: 



S 



d X\J—g 



R 



1 



- -d^d^ - - - V(4>) - U( X ) - (6(0) + 9(x)) g 



(111) 



This kind of action often appears in the models inspired by the string theory [Til ]. In such models, one scalar <j> may 
correspond to the dilaton and another scalar \ to modulus. We now consider the case that the derivative of \i ^^Xi 
is small or e is very small. Then we may neglect the kinetic term of x an d X could be regarded as an auxiliary field. 
Repeating the process (|92I94|) . we obtain the -F(G)-gravity coupled with the scalar field cf>: 



S 



d x\f—g 



^ - \d^d^ - V{cj>) - a(0)G + F(G) 



(112) 



The relation between scalar-Gauss-Bonnet gravity and modified Gauss-Bonnet gravity (or two parameterizations of 
the same theory) is discussed in this section. It is shown that cosm olog ical solutions obtained in one of such theories 
may be used (with different physical interpretation , compare with 27]) in another theory. It is often turns out that 
it is easier to work with specific parametrization of the same theory. Of course, only comparison with observational 
data may select the truly physical theory in correct parametrization. 



IV. CONCLUSION 



In this paper we have studied several aspects of (dilaton) gravity in the presence of string corrections up to third 
order in curvature. The second order term is Euler density of order two called, the Gauss-Bonnet term. The next-to- 
leading term contains higher order Euler density (-E3) plus a term of order three in curvature. The expression of E% 
is identically zero in space-time of dimension less than six; the term beyond the Euler density contributes to equation 
of motion even for a fixed field <fi. We have verified that the de-Sitter solution which exists in the case of Type II 
and Bosonic strings is an unstable node. It is shown that in the presence of a barotropic fluid (radiation/matter), 
inflationary solution exists in the high curvature regime for constant field. 

For a dynamically evolving field canonical in nature, there exists an interesting dark energy solution (|35|l char- 
acterized by H = ho/t, 4> = cf)Qlnt/ti for h > {H = ho/t s — t, <f> = 4>q ln(t s — t/ti) (when ho < 0). The three 
years WMAP data taken with the SNL survey [22] suggests that wde — — 1.06ltolol- W e have shown that choosing 
a range of parameter <pQ (which is amplified thanks to third order curvature term contribution) we can easily obtain 
the observed values of wde for phantom as well as for non-phantom dark energy. We have demonstrated, in detail, 
the stability of dark energy solution. For non-phantom energy, in the large ho limit, we presented analytical solution 
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which shows that one of the eigenvalues of the 3x3 perturbation matrix is real and negative where as the other two 
are purely imaginary, thereby, establishing the stability of solution l|35|) . We have verified numerically that stability 
holds for all smaller and generic values of kg in this case. The phantom dark energy solution corresponding to ho < 
turns out to be unstable. It is remarkable that string curvature corrections can account for late time acceleration and 
dark energy can be realized without the introduction of a field potential. 

It is shown how scalar-Gauss-Bonnet gravity may be reconstructed for any given cosmology. The corresponding 
scalar potentials for several dark energy cosmologies including quintessence, phantom, cosmological constant or oscil- 
latory regimes are explicitly found. This shows that having the realistic scale factor evolution, the principal possibility 
appears to present string-inspired gravity where such evolution is realized. It is explained how to transform scalar- 
Gauss-Bonnet gravity (even with account of third order curvature term) to modified Gauss-Bonnet gravity j2(| which 
seems to pass the Solar System tests. 

Different forms of modified gravity are attempted recently (for a review, see Q ) to describe dark energy universe; 
these models provide a qualitatively simple resolution of dark energy /coincidence problems and deserve further con- 
sideration. It is quite likely that time has come to reconsider the basics of General Relativity at the late universe in 
the search of realistic modified gravity/dark energy theory. 

We should also mention that in the present study we have tested the background model against observations. 
The study of perturbations in the scenario discussed here is quite complicated and challenging and in our opinion it 
deserves attention; we defer this investigation to our future work. 
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